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SUMMARY 


The  transfer  theory  is  applied  to  the  problem  of  atmospheric  diffu- 
sion of  momentum  and  heat  induced  by  line  and  point  sources  of  heat  on  the 
surface  of  the  earth.    In  order  that  the  validity  of  the  approximations 
of  the  boundary  layer  theory  be  realized,  the  thiokness  of  the  layer  in 
which  the  temperatures  and  velocities  differ  appreciably  from  the  values 
at  infinity  must  be  small  compared  to  the  height  of  the  region  of  influ- 
ence.   Variations  of  four  dependent  variables- -inflow  rate,  inflow  veloc- 
ities, maximum  velocities,  and  maximum  temperature--as  a  function  of 
elevation  for  reppesentative  heat  source  strengths  are  compared  for  the 
four  cases  investigated.    In  addition,  the  shape  of  the  region  of  influenoe 
and  the  factors  affecting  the  shape  are  investigated. 

Results  indicate  that  laminar  sources  tend  to  become  uniform  in 
cross-section  with  large  vertioal  velocities  and  temperatures,  thereby 
indicating  a  low  lateral  tendency  to  diffuse  bouyanoy  and  momentum.  The 
turbulent  sources  diverge  with  a  rapid  diffusion  of  bouyanoy  and  momentum 
depending  upon  the  mixing  length  proportionality  constant. 

Application  of  the  turbulent  point  source  results  to  a  finite  area 
source  has  been  presented  by  postulating  that  the  dependent  variables  of 
indraft  velocity,  maximum  temperature  and  maximum  vertioal  velocity  of  the 
finite  area  source  can  be  produced  by  a  point  souroe  of  equivalent  heat 
strength  functioning  at  a  specified  distanoe  below  the  finite  souroe. 
Data  concerning  the  behavior  of  the  aforementioned  variable!  induced  by  a 
finite  area  heat  souroe  do  not  exist.   Therefore,  experimental  investigations 
are  necessary  to  provide  guidanoe  for  further  analysis. 
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INTRODUCTION 


A  zone  of  convection  is  clearly  evident  when  a  fire  "burns  in  the 
free  atmosphere.    This  convection  zone  is  of  major  interest  in  understand- 
ing the  behavior  of  forest  fires.    The  transport  of  mass,  heat  and  momentum 
occurs  within  the  convection  column.    Knowledge  of  these  transport  processes 
will  provide  some  explanation  of  fire  behavior.    Quantitative  measures  of 
fire  behavior  must  stem  from  basic  explanations  of  behavior  which  are  con- 
sistent with  the  theories  of  heat  and  mass  transfer,  fluid  flow,  and  the 
chemistry  of  combustion  of  solid  fuels.    The  combustion  and  convection 
zones  are  inter-related.    The  oxygen  required  for  combustion  must  be  sup- 
plied by  the  convective  velocities.     In  turn  the  convection  zone  is  induced 
by  the  heat  liberated  by  the  combustion  process.    Despite  this  interaction, 
the  convection  zone  may  be  studied  by  utilizing  heat  transfer  and  fluid 
flow  theories.     If  a  quasi- steady  state  process  is  postulated,  the  mechan- 
ism governing  convective  flow  is  amenable  to  analysis.    The  present 
investigation  explains  the  free  convective  mechanism. 

Density  differences  which  are  associated  with  temperature  differ- 
ences due  to  thermal  expansion  cause  motion  of  a  fluid  which  is  called 
natural  flow.    The  natural  flow  Is  termed  free  convection  if  it  is  not 
subjected  to  constraining  boundaries.    Such  a  phenomena  exists  when  a  fire 
occurs  in  the  atmosphere  under  the  influence  of  negligible  wind  velocities. 
In  free  convection  the  velocity  and  temperature  distributions  are  inter- 
dependent, necessitating  the  simultaneous  solution  of  the  momentum  and 
energy  equations.    The  cases  investigated  here  are  for  the  condition  of 
zero  velocity  and  uniform  potential  temperature  in  the  ambient  region. 
The  effect  of  variable  potential  temperature  on  the  basic  system  of  equa- 
tions is  presented  for  the  laminar  line  source  by  an  approximate  method. 

The  earliest  study  concerning  the  effect  of  a  heated  source  was 
undertaken  by  Schmidt  (4>_5)— ^  in  1940.    He  presented  theoretical  investi- 
gations for  the  turbulent  line  and  point  sources.    In  addition,  he  obtained 
experimental  velocity  and  temperature  profiles  for  a  finite  electrically 
heated  source  55  mm  in  diameter  by  use  of  a  hot  wire  and  thermocouple 
traversing  probe.    Correlation  of  the  theoretical  and  experimental  results 
was  obtained. 

Rouse  et  al  (3)  conducted  experiments  to  determine  velocity  and 
temperature  profiles  for  line  and  point  sources.    The  two  dimensional 


1/    Underlined  numbers  in  parentheses  refer  to  Literature  Cited, 

Page  45. 
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investigation  was  confined  "between  parallel  walls  4  feet  high,  8  feet 
long  and  4  feet  apart  with  the  source  across  the  midsection  of  a  low 
platform  extending  the  length  of  the  walls.    The  axisymmetrical  case  was 
conducted  in  a  room  25  feet  in  diameter  with  an  over- all  height  of  11  feet, 
wherein  the  source  was  located  at  the  center  of  a  platform  8  feet  in  diam- 
eter o    Measurements  for  temperature  and  velocity  were  obtained  by  means  of 
a  copper- cons tan tan  thermocouple  and  a  vane  anemometer.    Recessed  gas 
burners  provided  rather  low  heat  source  output. 

Yih  (T>8)  presented  theoretical  investigations  for  laminar  free 
convection  due  to  line  and  point  sources  of  heat  for  small  variations  of 
pressure  and  temperature. 

The  present  investigation  presents  a  detailed  review  of  the  theo- 
retical analysis  of  laminar  and  turbulent  heat  sources  as  effecting  a  free 
convective  flow.    The  final  objective  is  to  establish  a  working  model  for 
the  convection-column  behavior  induced  by  a  fire  on  the  surface  of  the 
earth.    The  basic  analysis  presented  in  Appendix  A  has  been  formulated 
with  micrometeorology  in  mind  and  with  heat  source  strengths  representa- 
tive of  moderate  forest  fires.    Appendix  B  presents  the  investigation  of 
the  horizontal  pressure  gradient  since  this  gradient  is  taken  as  zero 
according  to  the  boundary  layer  theory. 

The  mathematical  investigations  are  essentially  of  a  simple  and 
idealized  nature  which  by  no  means  solves  the  complex  problems  of  the 
behavior  of  fire- convection  columns. 

RESULTS 

Since  the  solutions  obtained  here  do  not  account  for  variable 
properties,  the  magnitudes  involved  in  the  results  were  taken  at  an  average 
ambient  condition  for  air  at  atmospheric  pressure.    Thus  by  choosing  the 
properties  for  the  equations  and  selecting  reasonable  heat  rates,  a  com- 
parative study  of  the  variables  was  undertaken.    The  property  values  taken 
were 

Voo  =  1.80  x  10~4    sq.ft. /sec. 
poo  =  2.34  x  10-3  lb.sec.2/ft.-4- 

Cp  =  0.24    BTU/lb.  °R 
0oo  =  530  °R 

The  dependent  variables  obtained  from  the  theory  are  presented  in 
table  1.  The  mixing  length  proportionality  constant  was  taken  as  0.0417* 
which  was  experimentally  determined  by  Schmidt  (_5) .  This  value  was  taken 
in  order  to  effect  the  comparison. 
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Table  1  -  Dependent  variables  for  laminar  and 
turbulent  line  and  point  sources 


Line  Source  Point  Source 


Laminar 
^/L  =  l.l2l/oa|2>l/5Z3/5 

um  =  6.252/«,^3zy8 


=  0.67^-^/5Z-3/5 


^  =  3.00<2)J/5z2/5 


Turbulent 
'A 


=  0.97  Voo0'?>z 


■Jm=  7.40i/0o^1/3 


^m-^c  =  44.44  ^(2)2/3  z- 
&  --  0.30 Z 


Laminar 


<T=  1.0 


^  =  6.00  V<x>  Z 

=-0.502/oo  ^'^Z-'^ 


Um  =  0.802/oo<T2 


&  =  12.00 


Turbulent 
=  0.51  Voo#]/3Z5/3 

^Z^OT/oo^'^z-'/a 
um  =  38.50  Voo4>V*Z-V3 


^-^=35.70^=- 


&  =  0.36  Z 


Comparison  01  the  variables  is  presented  in  figures  1  to  5  for  two 
values  of  heat-source  strength.     These  rates  of  1.74  and  1740  BTU/sec.-ft. 
for  the  line  sources,  and  1.74  and  1740  BTU/sec .  for  the  point  sources 
were  taken  as  representative  values  of  energy  release  as  experienced  by 
actual  fires,    Only  one  heat  output  has  been  indicated  in  some  of  the  curves 
since  the  omitted  heat  output  yields  either  too  large  or  small  a  magnitude 
of  the  quantity  concerned. 
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Table  1  -  Dependent  variables  for  laminar  and 
turbulent  line  and  point  sources 


Line  Source  Point  Source 


Laminar 


er=  0.66 


%L--'l2W/5z3/5 
=-0.672/Oo//5z"2/5 

6.252/00 /^z'/s 


m 


0.67  2^L//5z-3/5 


^  =  3.00(2)J/5z2/5 


Turbulent 
'A 


0£/L  =  03lVoo<t>'*z 


u  m=  ^Oz/^'/s 


^m-#«,  =  44.44  ^^2/3z-' 
=  0.30  Z 


Laminar 


<7-=  1.0 


^  =  6.00  Voo  Z 

=  -0.502/oe  ^4  Z_l/2 


Um  =  0.802/00(^2 


<9m-^oo  =  0"  g 


^  =  12.00  (0-l/4  Z!/2 


Turbulent 

^  =  0.51  Veo4>X/*Zb/l 

um  =  38.50  v^^iz-^l 

em-e„  -  35.70  J^2/3z-% 

^  =  0.36  Z 


Comparison  of  the  variables  is  presented  in  figures  1  to  5  for  two 
values  of  heat-source  strength.    These  rates  of  1.74  and  1740  BTU/sec.-ft. 
for  the  line  sources }  and  1.74  and  1740  BTU/sec.  for  the  point  sources 
were  taken  as  representative  values  of  energy  release  as  experienced  by 
actual  fires.    Only  one  heat  output  has  been  indicated  in  some  of  the  curves 
since  the  omitted  heat  output  yields  either  too  large  or  small  a  magnitude 
of  the  quantity  concerned. 


9 


8.0 


• 

r~ 

UJ 

UJ 

u. 

to 

O 

z 

4.0 

o 

> 

UJ 

-I 

UJ 

M 

2.0 

0.0 
0.0 


7 

^— ALL  Q  BTU/SEC  / 
\  / 

-  \      /  7 

LAMINAR 

'  1 

 LINE  SOURCE 

 POINT  SOURCE 

\    /  / 
\  /  /x 

V  / 

^— Q/L  «  1.74 

BTU/SEC  FT 

/  / 
/  / 

/  / 

BTU/SEC  FT 

/  / 

^-0VL^=  1740 

/  s 

—  T"  I 

1 

1  i 

8.0 


4.0  8.0  12.0  16.0  20.0 

jjj  STREAM  FUNCTION  (FT3/SEC) 

 1  1  


24.0 


TURBULENT 


  LINE  SOURCE 

 POINT  SOURCE 


4.0  8.0  12.0  16.0  20.0 

STREAM  FUNCTION  (I03  FT3/SEC) 

Figure  1. --Theoretical  vertical  variation  of  stream 
function  for  laminar  and  turbulent  line 
and  point  sources. 


24.0 


10 


0.0  10.0  20.0  30.0  40.0  50.0  60.0 

<?m-S«>  MAXIMUM  TEMPERATURE  RISE  (°F) 


0 


0 


0 


1             1  1 
TURBULENT 

1 

  LINE 

 POINT 

I 

SOURCE 
SOURCE 

^ —  q/l  = 

1740  BTU/SEC  FT 

z-Q  =  1740  BTU/SEC 

V  / 

i 

0.0  50  100  150  2  00  250  3  00 

0m-5oo  MAXIMUM  TEMPERATURE  RISE  (<>F) 


Figure  4 • --Theoretical  maximum  temperature  variation 
for  laminar  and  turbulent  line  and  point 
sources . 


13 


$  ACTION  ZONE  (10"'  FEET) 


6.0 


&  ACTION  Z0NE(I03  FEET) 


Figure  5 . --Theoretical  variation  of  action  zone  for 

laminar  and  turbulent  line  and  point  sources . 


DISCUSSION 


Although  these  investigations  of  the  laminar  hypothesis  and  the 
Prandtl  mixing- length  hypothesis  are  idealized,  they  are  readily  amenable 
to  mathematical  solution.    The  approach  has  been  to  initially  retain  math- 
ematical simplicity,  thereby  avoiding  complex  differential  equations. 
After  the  idealized  investigations  provide  a  feeling  for  the  problem, 
modification  of  the  system  is  necessary  to  conform  with  the  physical  sys- 
tem, provided  that  data  are  available.    Application  of  the  results  to 
actual  fires  must  be  done  by  keeping  in  mind  that  the  theoretical 
development  has  been  formulated  according  to  the  following  conditions: 

1.  Steady  state 

2.  Constant  properties 

3 .  Omiss  ion  oi  c  ombustion  zone 

4.  Omission  of  radiation  heat  transfer 

5.  Dry  air  (perfect  gas) 

6.  Zero  velocity  and  uniform  potential  temperature  in 
ambient  region 

Certain  of  these  conditions  present  serious  limitations  when  con- 
sidering application  of  the  theory  to  the  real  problem  of  fire  behavior. 
The  combustion  process  and  the  transient  nature  of  the  actual  problem  are 
important.    Combustion  of  a  forest  fuel  may  occur  over  a  range  of  combus- 
tion efficiency  which  is  believed  to  be  controlled  by  chemical  reactivity 
and  the  diffusion  of  oxygen.    The  slower  process  will  determine  the  combus- 
tion rate.    A  change  in  the  combustion  process  gives  rise  to  the  transient 
problem.    As  a  consequence  the  rate  of  combustion  attains  a  new  equilibrium 
level. 

The  analysis  has  considered  the  vertical  pressure  gradient  in  terms 
of  the  density  of  the  exterior  region.     For  the  energy  equation  the  verti- 
cal pressure  gradient  has  been  accounted  for  by  the  potential  temperature. 
The  potential  temperature  is  then  utilized  to  eliminate  the  density  dif- 
ference in  the  momentum  equation.    This  potential  temperature  difference 
generates  the  velocities .    What  is  implied  by  variable  properties  is  the 
influence  of  the  variation  of  the  absolute  viscosity  and  density  from  the 
exterior  region  to  the  center  of  the  convection  column.    The  analysis  of 
this  condition,  that  is,  one  considering  compressibility,  remains  mathemat- 
ically intractable  at  the  present  time. 

In  actual  combustion  of  a  solid  fuel,  water  vapor  arises  from: 
the  moisture  driven  from  the  fuel,  the  moisture  resulting  from  the  oxida- 
tion of  the  fuel,  and  the  moisture  existing  in  the  ambient  air.    The  con- 
centration of  water  vapor  in  the  convection  column  could  be  determined  if 
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the  oxidation  reaction  for  the  solid  fuel  is  known.    In  the  convection 
column  the  water  vapor  provides  a  region  of  nonuniformly  distributed  heat 
sources  which  induce  a  velocity  field  similar  to  a  thunderstorm  cell  in 
the  cumulus  stage  of  development. 

Considering  conservation  of  mass  for  constant  properties ,  the  rate 
of  change  of  velocity  into  the  column  with  elevation  is-- 


2L& 


o  um  d(  JO 


d(um<?) 
dz 


line  source, 


2^ 


''_y.iL  Ai. 


d(Um^2) 
dz 


point  source, 


Considering  conservation  of  energy  which  is  produced  "by  the  heat 
source  for  constant  properties  yields- - 


Q_ 
L 


u    a  —ea 


(f) 


line  source, 


*'  y  u  (e-e<»\d 


&rr\~&<»)  point  source, 


The  foregoing  relations  are  valid  for  finite  sources  provided  that 
velocity  and  temperatures  profiles  are  universal  in  character.    If  the 
equations  for  the  line  and  point  sources  are  utilized  in  the  conservation 
of  mass  relations,  the  integrals  are  constants.    For  the  conservation  of 
energy  relations,  substitution  of  the  line  and  point  source  equations  for 
the  turbulent  cases  yields  constants  which  are  dependent  on  the  mixing 
length  proportionality  factor. 

Comparison  of  these  cases  of  laminar  and  turbulent  flow  indicates 
diverse  characteristic  behavior.    Figure  1  shows  the  vertical  variation 
of  the  stream  function  from  which  the  indraft  rate  may  be  obtained.  If 
the  velocity  and  temperature  distributions  are  universal,  the  mass  and 
energy  balance  relations  indicate  that  the  rate  of  change  of  the  inverse 
of  the  maximum  temperature  difference  with  elevation  controls  the  indraft 
velocity i    For  instance,  in  the  line  source  case,  if  the  maximum  tempera- 
ture difference  decays  as  the  inverse  of  the  elevation,  then  the  indraft 
velocity  is  constant.    From  these  relations  several  combinations  of  the 
dependent  variables  are  possible.    Only  experimental  evidence  concerning 

universal  nature  of  the  velocity  and  temperature  profiles  and  the 
behavior  of  the  independent  variables  can  provide  the  final  answer. 
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Comparison  of  the  entraining  velocities  is  presented  in  figure  2, 
which  shows  a  decrease  in  velocity  with  an  increase  in  elevation,  except 
for  the  turbulent  line  source,  which  is  constant.    Laminar  velocities  are 
smaller  in  magnitude  than  the  turbulent  velocities  for  the  same  source 
strength.    Turbulent  velocities  are  more  sensitive  to  the  source  strength 
than  laminar  velocities. 

Maximum  vertical  velocities  which  occur  at  the  perpendicular  axis 
of  the  sources  are  presented  in  figure  3°  The  higher  heat  output  for  the 
laminar  case  has  not  been  shown  since  the  magnitudes  are  above  the  veloc- 
ity of  sound.  These  results  indicate  greater  velocities  for  the  laminar 
case  than  the  turbulent  for  a  constant  source  strength.  In  addition,  the 
laminar  source  velocities  are  more  sensitive  to  source  strength. 

Axis  potential  temperature  variations  are  pre;  .ented  in  figure  4« 
Results  for  the  larger  source  output  for  the  laminar  cases  have  been 
omitted  since  the  temperatures  are  excessive.    In  contrast,  the  smaller 
heat  output  for  the  turbulent  cases  has  been  omitted  due  to  the  small 
resulting  temperatures.    In  all  cases  the  potential  temperatures  decrease 
with  height  which  is  in  agreement  with  the  physical  system. 

Figure  5  shows  the  shape  of  the  zone  of  influence  due  to  the  sources 
investigated.    The  laminar  sources  indicate  a  very  narrow  spread  becoming 
more  concentrated  as  the  source  output  is  increased.     In  direct  contrast 
the  turbulent  cases  expand  very  rapidly  in  a  linear  behavior,  and  are  in- 
dependent of  source  strength  on  the  basis  that  the  mixing  length  propor- 
tionality factor    c    is  a  true  constant .    Since  this  factor  does  appear  in 
the  majority  of  the  variables  of  the  turbulent  cases,  there  remains  the 
possibility  that  this  proportionality  factor  may  depend  on  other  parameters. 
Only  further  investigation  of  the  atmospheric  diffusive  mechanism  can 
determine  this  behavior. 

Recent  inquiries  have  been  directed  to  the  effect  of  lapse  rate  on 
the  convection  column  which  would  thereby  affect  the  heat  source.    In  Ap- 
pendix G  an  approximate  investigation  has  been  undertaken  for  the  effect 
of  nonuniform  potential  temperature  for  the  laminar  line  source.  This 
analysis  indicates  that  the  gradients  of  potential  temperature  must  be  very 
large  to  be  contributing  factors  for  heat  source  outputs  considered  in  this 
report.    This  conclusion  was  also  obtained  by  Byram  (l)  in  his  thermody- 
namic analysis.    For  small  heat  source  strengths  the  lapse  rate  is  undoubt- 
edly of  importance  as,  for  instance,  in  whirlwinds  or  cumulus  cloud  develop- 
ment.   The  foregoing  does  not  imply  that  the  lapse  rate  is  unimportant  in 
the  problem  of  atmospheric  fire  behavior  since  the  eddy  viscosity  and 
conductivity  are  greatly  affected  at  least  near  the  surface. 
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CONCLUSIONS 


The  present  investigation  shows  that  the  laminar  sources  yield 
narrow  convection  columns  which  consequently  require  large  vertical  veloc- 
ities and  temperatures,  extending  the  column  to  higher  elevations.  Since 
the  turbulent  sources  spread  in  a  horizontal  direction,  the  vertical 
velocities  and  temperatures  are  smaller  in  magnitude,  and  the  column  does 
not  attain  comparable  heights. 

From  the  results  which  indicate  the  magnitudes  of  the  variables, 
it  may  be  concluded  that  the  laminar  process  is  not  compatible  with  the 
physical  system.    Therefore,  the  results  of  the  turbulent  investigation 
should  be  used  for  practical  application. 

A  source  of  heat  induces  a  velocity  field  wherein  momentum  and 
heat  must  be  diffused.    For  the  condition  of  negligible  free  wind  speeds, 
the  mode  of  diffusion  plays  an  important  role  in  the  behavior  of  the 
variables.    This  may  be  understood  in  that  the  region  of  influence  of  the 
heat  source  controls  the  velocity  and  temperature  profiles  and,  by  virtue 
of  conservation  of  mass  and  energy,  controls  the  indraft  requirement. 
Characteristics  of  the  heat  and  momentum  diffusion  must  be  determined  to 
ultimately  understand  the  problem. 

For  large  heat  outputs  the  potential  temperature  gradients  do  not 
influence  the  results  of  the  basic  theory  for  an  established  transfer 
process,  but  are  important  in  effecting  the  diffusion  mechanism. 


RECOMMENDATION 

An  experimental  investigation  to  determine  the  effect  of  radius 
and  heat  rate  on  the  indraft  velocities  and  velocity  and  temperature 
distributions  in  the  convection  zone  is  recommended. 2/ 

This  experiment  should  be  conducted  in  the  laboratory  so  that  the 
eddy  viscosity  and  eddy  conductivity  would  not  be  affected  by  the  lapse 
conditions.    Then  the  condition  of  linear  mixing  length  with  height  could 
be  fulfilled  so  that  the  mixing  length  proportionality  constant    c  would 
be  a  true  constant.    The  system  should  be  designed  for  actual  flame  sources 
whereby  the  heat  production  rate  could  be  controlled  for  a  specific  range 


2/    Application  of  the  turbulent  point  source  to  the  case  of  a 
finite  area  source  of  heat  is  presented  in  Appendix  D. 
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of  values.    The  sizes  of  the  flame  front  would  be  determined  by  the  maxi- 
mum energy  release  which  is  feasible  such  that  heat  rates  of  sufficient 
magnitude  produce  measurable  velocities. 

Also  the  system  should  be  flexible  so  that  electrical  heaters  could 
replace  the  flame  sources  in  order  to  find  the  effect  of  sources  with  and 
without  combustion.    Shielded  thermocouples  would  be  designed  for  tempera- 
ture traverses,  and  the  heated- thermocouple  type  anemometer  and  hot  wire 
anemometer  for  velocity  measurements.    A  traversing  mechanism  would  allow 
the  two  measuring  circuits  to  be  located  in  the  vertical  and  radial 
directions. 

Damping  screens  located  at  ample  distances  surrounding  the  source 
would  provide  the  necessary  damping  so  that  any  exterior  disturbances  would 
decay  thereby  not  influencing  the  convection  column.    This  would  also  allow 
the  indraft  makeup  to  feed  the  convection  column.    The  height  of  the  damp- 
ing section  would  be  determined  by  the  growth  of  the  convection  zone. 
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APPENDIX  A 


BASIC  ANALYSIS 


Laminar  Line  Source 

The  continuity,  momentum,  and  energy  equations  for  the  two- 
dimensional  flow  due  to  a  line  source  for  a  fluid  of  constant  Prandtl 
Number,  specific  heat,  and  viscosity  are 

fr+fr°  cd 

du       du       dp      d  u  /  x 

pgJCpufI+MJCpv|I-u^-=  Jk|i  (3) 
v  dz  v  dy      dz  dy2 

dp 

wherein  the  condition  that       is  of  order  zero  has  been  obtained  from  a 
non-dimensional  order  of  magnitude  investigation  of  the  Navier-Stokes 
equations.    Applying  equation  (2)  in  the  region  exterior  to  the  region 
of  influence  there  results: 

dz     Hc°q  dz 

since  p  =      .       The  foregoing  equations  become 
du  dv 


dz  dy 


=  0  (4) 


du         du      ,        x  dZu 


dJ  dT      dp      ,  d2T 


pgJcDu- —  +  pgJcnv-  u^-^Jk-r— 5-  (6) 

p  dz       *    P  dy       dz  dy2 


From  the  perfect  gas  law  at  constant  pressure 
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which  gives 

Using  the  definition  of  potential  temperature 
I    de    i  fdT        i  dp\ 


and 


0  dz    T  \dz     Jcppg  dz 


\   de     \  dl  i  <?2e?     i  d2T 


e  dy     T  dy  e  dy2  '  T  dyz 

the  foregoing  equations  are 

du     dv    ^  ,  x 

+  -7^7  =  0  (7) 


c?z  <?y 


u —  v —  =  —s   (9) 

dz  dy  cr  dyd 

with  the  boundary  conditions 

y-0     v-0     -^  =  0   -0 

y— *i  00     u  =  0      0 -^co  =  0 

The  equations  have  been  derived  on  the  basis  that  the  density  difference 
is  the  important  factor  compared  to  the  magnitude  of  the  density  itself. 
Thus  the  variation  of  the  density  and  absolute  viscosity  with  temperature 
has  not  been  considered. 


Solution  of  equations  (8)  and  (9)  proceeds  as  follows 
Q 


'/5 


\-cpea>pO0-u^>J  z2/5 


(10) 
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(11) 


so  that 


U  =  —r-  -  Zf 


9  f5z'/5^- 

3  /    z  d?7 


(12) 


Q 


(13) 


e  -em  = 


where  u     and  v    as  so  chosen  satisfy  the  continuity  equation  (7). 


(M) 


In  addition,  the  condition  that  the  heat  output  of  the  source  is 
equal  to  the  heat  flow  at  any  cross-section  of  the  region  of  influence 
yields 


Q 


p^q  Cpu(e-<9oo)  dy 


(15) 


With  the  foregoing  transformations,  equations  (8),  (9)  and  (15) 
are  transformed  into  a  system  of  ordinary  differential  equations  as 


5^  +  3fd2f 


df 

d?73  '      d??2  Vd?7 


+  5h  =  0 


(16) 


5  d2h  d(fh)  _  n 
3<T  d?72      d/7     '  U 


(IT) 


+■00  At 

d77  =  ! 

-00  d?7 


(18) 


with  the  corresponding  boundary  conditions 
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f(0)  =  f"(0)  -  h'(0)  =  0 


f'(ioo)  =  h(t»)  =  o 

For  a  Prandtl  Number  of  2/3  the  functions  satisfying  the  foregoing 

equations  and  boundary  conditions  are 

f(77)^  tanh  ^-7]  (19) 


h(77)--^sech2  ^f-r,  (20) 

Turbulent  Line  Source 

The  equations  for  this  case  are  written  by  utilization  of  the 
Prandtl  mixing  length  hypothesis  wherein  the  eddy  diffusivity  for  heat  and 
momentum  are  considered  equivalent  and  as  JJ2^L         thereby  yielding 

^  +  £■=0  (21) 
dz  dy 

\2d 


du  du  \ 
dy  dy  ) 


u|u^v^  =  _^V  dy  d_y±  fsp^ 
dz      dy  dy  *\&oo 

du  die-eooT 

u  die-ej  |  yd{e-eJ  _    dV  dy     dy     i  (23) 
dz  dy  dy 

where  4  =  cz 

It  should  be  pointed  out  that  the  dependent  variables  used  herein 
represent  mean  quantities. 

By  use  of 

77=1-  (24) 
=  ^)3zf(7)  (25) 


so  that 


dijt 

u  =  —r-  -  va 


Q 


df 


0  -©o 


Q 


?/3 


zH  h(?7) 


yields  the  transformed  equations  and  condition  as 


2  \2 

df 


:2dv^__f^!i  _0 

d772  d  772 


d  f  dh 

c2jdVd77?_dW     d(hf)  = 
d  77  d  77 


with  boundary  conditions 

f(0)  =  f"(0)  =  h'(0)  =  0 

f'(ico)  =  h(lco)  =  0 

?  A 

Letting      77  =  (2  (T )  ^77,         and       h  = 
and  (30)  become 


H 


(2c2f/3 


dT/2    dT/3         d  77j^ 


2? 


67]  d77 

The  functions  f    and  H    have  "been  obtained  by  Schmidt  (4)  who 
utilized  a  series  solution  approach  wherein  the  first  set  for  f    and  H 
applied  near  the  region  y  =  0    and  the  second  set  near  the  edge  of  the 
region  affected  by  the  source. 

The  set  near  the  central  axis  (y  =  0)  is 

f(77,)  =  I.I39877J-  0.3771 77^/2 +  0.0633 77,4- 0.0051 771"/2+  

(34) 


H(77)  =  I  -  l.0746773/2  +  0.6I0877,3  -  0.1154  77^ 


This  approach  is  utilized  since  the  foregoing  series  do  not 
converge  adequately  for  large  values  of  77  ( . 

The  variables  become 

U-CpS.fl.ig/     (2c2)'/3  dT?, 


e-e~- 


lie*  (       Q       V3  H 


9    VLCp^ooPoo^/  {2c2fz 

Laminar  Foint  Source 

For  this  case  the  equations  for  the  ax i symmetrical  flow  are 

^>+^l=0  (36) 
az  ay 
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dz  y        <?y  y 


r  ^(g-gJi 

tfz  cty         cry  <?y 

By  use  of 


Q     V/4  V 


so  that 


V-—L<X- -^(—S—J*     (2±_  «) 

yields  the  transformed  equations  and  conditions  as 


77-^- 

d(hf)  _  _  j_  dV  d?7 
d  77    ~    4<7  d?7 


2? 

I 


df 


with  boundary  conditions 

f  (0)  =  f"(0)  =  h'(0)  =  0 

f'(i<»)  =  h(i«0  =o 

For  a  Prandtl  Number  of  unity  the  functions  satisfying  the  fore- 
going equations  and  boundary  conditions  are 

fw=l('  ! — r) 

V  I   +   r-   V 

Turbulent  Point  Source 


are 


The  axisymmetrical  equations  utilizing  the  Prandtl  mixing  length 


^  +  4^  =  0  (49) 
dz  dy 


(.2.  du  du\ 
du       du     i    dv  y  dy  dy    ,  fe-ea>\ 
z       dy     y  cty  \  g»  / 


(50) 


r.2         g-gJ  i 

• 1  y^y-^ — J 


c?(g-gj  ,  „  die-ej  _    i        y  dy  dy 

— ~3  hV  3  T   151J 

dz  dy  y  dy 


2* 


where  J|  -  cz 

By  use  of 

y 

— 


>p&coP< 


so  that 


<7> 


2/0 


\*  |_  df  z../3 


v'    y   tfz  "  Mcpft^z/J 


e  -ft, 


2/oofta 


Q       \2/3  „ 


9  ^Cpftcft,^ 


■5/3  h(77) 


yield  the  transformed  equations  and  condition  as 

\2- 


didf 


7?  d?7 
d77 


d77 


__J_±(dA  _  5_f  d 


J_  dl_\ 
77  d??/ 


3  77  \d77/     3         d 77 


f  (L^~) 

dh 

d 

[  d77 

d>7. 

5  d(fh) 

d/7 

3  dv 

h^dr?--\ 
-00    d  77 


with  the  boundary  conditions 


f(o)  =  f"(o)  =  h'(o)  =  0 

f'U-)  =  g(i-)  =  0 


H 


Letting  77  =  (3C  and  h=  (3C2)/3      equations  (57)  and  (58) 

are  the  same  as  those  derived  by  Schmidt  (_5) .  Again  using  the  series 
development  as  performed  by  Schmidt  there  obtains  for  the  axis  region 


f(77()  =  0.5405  7?f  -  0.275077^2  +  0.0642  -  0.008677'3/2  +  ... 
H(77)  =  I- 1.2480  >7^2  +  0.6802  773  -0.2220/79/2  +  0.0500  n6 


yielding 


CD 


v'/3 

Q 


1J)  z% 


'  Cp  @oo  Poo  Uoo 

U-_VJ  S  f  _!        z-h  _LJli 

Q     f3    i     z-'/3  dt 


'I 


(60) 


(61) 


m 


APPENDIX  B 


INVESTIGATION  OF  HORIZONTAL  PRESSURE  GRADIENT 


Since  the  analysis  presented  herein  is  based  upon  the  approxima- 
tions of  the  boundary  layer  theory  wherein  the  pressure  gradient  in  the 
y    direction  is  taken  of  zero  order,  an  investigation  of  the  validity  of 
the  approximation  is  presented  for  the  line  sources. 

Taking  the  momentum  equation  in  the  y    direction  as 

U7I+]'W~'i'~to~  (Laminar)  (62) 

Transforming  the  equations  from  z,y    to  77  there  results 


n      n    -  ft>^% 
Pm  4,  

0     m  % 


'  (9fS  .2  5V5"  „  .     h  5V5"  „ 

y   Q      sech — - — 77  tanh— - — 77 


^  secn<  iVF 375  ^  sech2  5VT  ^  tQnh  ^vyv  (63) 


2lV5    .    .    5V5"  „     5   „  ea^h2  5V5 

tanh      — 77-— 77  sech  — ^—V)dq 


z2    Jo  V  125 


Using  the  largest  value  of  the  heat  output  function  used  in  the 
previous  results  as  </)  -  I015  and  the  property  values  as  previously  used  of 
pM  =  2.34xl0~3  and  V-  =  1.80  x|Q"4  there  results 


7.6  xlO"5  t   .  7.6x10"" 
Z' 


r.t>  x  u      r      r.bxiu      .  t,  ,v 

P^"Pm"  3   !2 


Thus  even  for  small  values  of  z    the  pressure  drop  is  of  zero 

order. 

For  the  turbulent  line  source  the  result  given  by  Schmidt  (4)  is 
Pcf-Pm  =  0.39(2C2)2/3^^2/3  (65) 

3X 


Using  the  previous  values  of  the  quantities  and  taking  c=  0.0417 
there  results 

-  pm  =  0. 00675  lb./sq.  ft.  -  zero  order 
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APPENDIX  C 


EFFECT  OF  NONUNIFORM  POTENTIAL  TEMPERATURE 


To  obtain  an  idea  as  to  the  effect  of  vertical  potential  temper- 
ature gradients  in  the  atmosphere  upon  a  heat  source,  an  approximate 
approach  was  undertaken.    The  analysis  shall  be  presented  for  the  laminar 
line  source. 

Taking  the  continuity,  momentum,  and  energy  equations  as  previously 

derived 


and  use  of  the  Mises  transformation. 


d\j  v      **>  dz 

which  satisfies  the  continuity  equation  (66)  if  p^  -  constant.  Equations 
(67)  and  (68)  become  with  the  condition  that  fi Q  - U^p^  3/ 


du 

...  &£djl  .  9 
dz     00    dili  u 


0{zj)  , 


(69) 


cte  .  ^  d{u%)  ,  . 


_3/    This  relation  adequately  describes  the  variation  of  viscosity 
of  air  with  temperature  as  outlined  by  Chapman  and  Rubesin  (2). 


3> 


By  assuming  velocity  and  temperature  profiles  which  satisfy  the 

"boundary  conditions  of  the  problem  as 

.£ 

u{zj)-  um(z)  e  d*  (71) 

where  u(z,oo)=0         and  u'(z,0)  =0 

and 

<9(Z,O)-0(Z,C»)  "e  {12) 

where  e  (z,oo)  .=  6>(z,co)  <9'(z,0)=0 

Upon  substitution  of  these  distributions  into  equations  (69)  and 
(70)  at  the  condition  tp-0     by  reason  that  the  central  axis  is  the  region 
wherein  the  greatest  variations  in  velocity  and  temperature  occur,  the 
final  equations  are 

dz      V    U  e(z,a>)U  l™ 

dB.dG.     2    U  n  , 

-  "  —  "T2  0  (  74) 


and 


dz     dz       C"  d' 


from  the  condition  that 


_Q 
L 


/-+00 

^gcpu[^(z,^)-e(z,c»)]  dy 


where 


0  = 


B 


6>(Z,0)  -g(z,oo) 

e{  0,00)2/! 

g 

^(z.co) 
^(0,00)  vt 


<t>-- 


Q 


Lcp  ^(0,00)^14 


By  eliminating  0    and  B   from  the  three  equations  there  results  a 
second  order  nonlinear  differential  equation  for  the  unknown  function  D 


dz  dz^  Vdz/ 
D  d2!"?) 


9_  J.  _ 


V7F   2/i  B 


7V7T 


D4^dDj  d(5) 


dz 


+  V77  D 


77"  D 


dz  dz2 
7  d\^ 


V77  D' 


5  d' 


:D  d(5) 


dz^  dz 


-  677-D 


6  dD 


dz 


li) 


dz 


dz 


dz' 


+  V77"  77"  D 


0) 


dz 


(76) 


with  the  boundary  conditions 

Z-0     D-0     since    ll(0,^)  =  0     and  0(o>)=(O,oo) 

z  —  00     D  — 00   since  u(oo,^)  =  um(oo)     and    e{co,i//)  -  0(co,O) 

As  a  preliminary  verification  of  equation  (76)  a  check  for  the 
constant  potential  temperature  case  can  "be  obtained.    For  this  case 

B  =  Const  =g/i/<2>  Thus  equation  (76)  is 


n3  dD  d2D  .  on2/dD\3       4  . 
~dz~  ~&z}  Vdz/  ~  ~7wr 


~frT 


which  has  the  solution  satisfying  the  boundary  conditions 


3? 


thereby  yielding 


(T7) 


These  results  are  in  the  same  form  as  obtained  by  the  rigorous 
solution  for  the  laminar  line  source  as  given  in  table  1.    The  constants 
do  not  compare  favorably. 

The  nonlinear  equation  (76) ,  which  does  not  have  a  formal  solution 
available,  could  be  solved  by  a  numerical  method  or  series  development  of 
D    in  terms  of  z    for  the  non-dimensional  parameter  B/^-    This  approach 
was  not  deemed  necessary  since  an  investigation  of  the  order  of  magnitude 
of  the  factor  B/^  discloses  a  magnitude  of 

1.73  x  IP'3  g(z,oo) 

Q  0(O,oo) 
L 

For  the  heat  outputs  considered  in  this  report,  gradients  of  poten- 
tial temperature  would  necessarily  be  very  large  such  that  the  terms  of  the 
right-hand  side  of  equation  (76)  become  comparable  to  the  terms  of  the 
left-hand  side.    Thus,  this  order  of  magnitude  argument  indicates  that  the 
potential  temperature  gradient  is  a  higher  order  effect.    This  does  not 
imply  that  the  lapse  rate  is  unimportant  in  the  problem  of  atmospheric 
fire  behavior. 


APPENDIX  D 


APPLICATION  OF  FINITE  SOURCES 


Investigation  of  the  finite  source  dependent  variables  from  an 
exact  mathematical  approach  presents  a  system  of  extreme  complexity.  In 
the  point  source  investigation,  the  boundary  layer  assumptions  provide  a 
single  momentum  equation  which  cannot  he  realized  in  the  finite  source 
case.    To  obtain  an  order  of  magnitude  of  the  variables  of  interest,  the 
point  source  analysis  has  been  applied  to  the  finite  area  source  by  making 
definite  assumptions. 

The  postulate  is  made  that  the  finite  heat  source  of  total  heat 
strength  of  Q   BTU/sec.  is  equivalent  to  the  point  source  of  identical  heat 
rate  functioning  at  a  fictitious  distance  z0  below  the  actual  finite  source 
wherein  the  equations  for  v^y,  ^-(9^      and  um    are  applicable.    The  foregoing 
statement  implies  that  there  exists  an  equivalent  point  source  which  yields 
the  identical  behavior  of  the  dependent  variables  as  the  actual  finite 
source.    The  vertical  distance  is  modified  to  make  the  equality.    Such  a 
condition  could  be  reasonable  for  the  region  exterior  to  the  immediate 
vicinity  of  the  finite  area  source.    Furthermore,  if  these  results  are  to 
be  applied  to  actual  fire  behavior,  the  reader  must  realize  that  the  com- 
bustion zone  and  flame  front  have  not  been  considered  in  the  free  convec- 
tion investigation.     In  addition,  actual  fire  behavior  is  a  transient 
phenomenon  from  build-up  to  decay  which  must  be  the  subject  for  final 
consideration. 


Since  data  on  finite  sources  of  heat  in  free  convection  are  virtu- 
ally nonexistent,  the  determination  of  the  parameters  which  influence  z0 
must  remain  simple  at  this  stage.    The  quantity  is  obtained  from  the 
boundary  layer  growth  of  the  point  source  which  is 

R_ 


z0=  o,,/.  (78) 


thereby  indicating  that  the  fictitious  distance  is  a  function  of  the 
radius  of  the  finite  source  and  the  mixing  length  proportionality  constant. 
It  is  fully  realized  that  other  parameters  do  influence  this  determination, 
Until  experimental  evidence  is  obtained,  no  refinements  are  considered 
warranted. 

By  the  reasoning  outlined  the  dependent  variables  become 
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=  (3C2)2/3  (81) 

2 

Substitution  for  z0      the  condition  that  Q=77"R  q  and  using 
the  property  values  previously  indicated  there  results 
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1.63  q'/3Rl/5 
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The  foregoing  equations  relate  the  dependent  variables  to  the  heat 
rate  per  unit  area,  size  of  heat  source,  and  mixing  length  proportionality 
constant .    At  this  stage  the  indraft  velocity  has  not  been  corrected  for 
surface  shear.    Investigation  of  the  mixing  length  proportionality  constant 
nas  been  conducted  by  Sheppard  (6).    This  study  indicates  that  the  mixing 
length  is  approximately  linear  with  height  only  for  the  neutral  case 
(adiabatic  lapse  rate).    Further  research  concerning  the  turbulent  mixing 
in  the  atmosphere  or  the  effect  of  turbulence  on  controlled  fires  need  be 
undertaken „ 

The  relations  derived  herein  were  applied  for  the  values  of 
q  =  20  BTU/see^ft„2  and  a  radius  of  100  ft.  for  values  of  C    equal  to  0.50 
and  0.30.    The  heat  rate  £aken  corresponds  to  a  forest  fuel  rate  of 
3  tons/acre-min.    The  results  indicate  that  the  indraft  velocities  are  not 
?ery  sensitive  to  the  mixing  length  proportionality  constant.    The  tempera- 
ture distribution  given  by  the  relation  is  for  a  convection  column  not 
accounting  for  the  flame  region.    Temperatures  in  the  flame  region  are 
greater  in  magnitude  than  the  values  given  by  the  relation.  Temperature 
distributions  in  flames  is  another  experimental  study  of  considerable 
interest.    The  vertical  velocities  are  influenced  considerably  by  the  mix- 
ing length  proportionality  constant.    These  results  have  been  plotted  in 
figures  6  and  7  respectively. 


A  model  for  the  convective  column  dependent  variables  generated 
by  a  finite  area  steady  state  heat  source  has  been  established  by  utiliz- 
ing a  basic  theoretical  approach.    In  order  to  obtain  the  reliability  and 
apply  refinements  to  the  model  advanced  herein,  an  experimental  investiga- 
tion must  be  conducted  with  the  mixing  length  proportionality  constant 
invariant  but  with  the  variation  of  the  heat  rate  and  size  of  the  flame 
front o    Instrumentation  is  required  to  obtain  the  dependent  variables. 
In  this  manner  the  logical  foundation  of  a  complete  model  for  ultimate 
fire  behavior  may  become  realized.    Continued  progress  from  the  theoretical 
approach  must  be  aided  and  guided  JJy  controlled  experiments  with  the  In- 
fluence of  the  heat  rate  and  size  and  source  providing  the  basic  foundation. 
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NOMENCLATURE 


The  following  nomenclature  is  used  in  the  paper: 

B  =  Function  of  equation  (74-) 

c  =  Mixing  length  proportionality  constant 

cp  =  Specific  heat  at  constant  pressure,  BTU/lb.  °F 

d  =  Function  of  equations  (71)  and  (72),  ft.2/sec. 

D  =  Non-dimensional  function  of  equations  (73) ,  (74)  and  (75) 

f(^)  =  Function  of  the  independent  variable 

g  =  Gravitational  constant,  ft ./sec. 2 

h (77)  =  Function  of  the  independent  variable 

H  ( 77( )  =  Function  of  the  independent  variable 

J  =  Mechanical  heat  equivalent,  ft.lb./BTU 

k  =  Thermal  conductivity,  BTU/sec- sq.ft.  °F/ft. 

i  =  Mixing  length,  ft. 

L  =  Line  source  length,  ft. 

p  =  Pressure  in  source  zone,  psf 

Pa.  =  Atmospheric  pressure,  psf 

Q  =  Heat  output  of  source,  BTU/sec. 

q  =  Heat  output  per  unit  area,  BTU/sec. -ft 

R  =  Radius  of  finite  source,  ft. 

T  =  Temperature  in  source  zone,  °R 

Too  =  Ambient  temperature,  °R 

u  =  Velocity  component  in  vertical  direction,  fps 

Um  =  Vertical  velocity  at  source  axis,  fps 
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U  =  Function  of  equations  (73)  and  (74),     l /ft. 

v  =  Velocity  component  perpendicular  to  source  axis,  fps 

y  =  Distance  normal  to  source  axis,  ft. 

z  =  Distance  vertical  from  point  source,  ft. 

Z1  =  Distance  vertical  from  finite  source,  ft. 

Zq  =  Fictitious  vertical  distance  for  finite  source,  ft. 

Tj  =  Independent  variable 

=  Turbulent  independent  variable 

6{z,fi)  =  Potential  temperature  for  approximate  method  appearing  in 
equation  (69),  °R 

(z,oo)  =  Atmospheric  potential  temperature  for  approximate  method,  °R 

(0,oo)  =  Surface  potential  temperature  for  approximate  method,  °R 

0  =  Function  of  equations  (73),  (74)  and  (75) 

/J  =  Absolute  viscosity,  lb. sec. /ft. 2 

2/oo  =  Ambient  kinematic  viscosity,  ft.2/sec. 

p  =  Mass  density  of  source  zone,  lb.sec.^/ft .4 

p^  =  Ambient  mass  density,  lb.sec.2/ft .4 

0~  =  Prandtl  Number,/^  Cp /k  non-dimensional 

0  =  Heat  output  function 

i)j  =  Stream  function 

0  =  Potential i temperature  in  source  zone,  °R 

e  =  Ambient  potential  temperature,  °R 

00 

0m  =  Maximum  potential  temperature  in  source  zone,  °R 

&  =  Action  zone,  ft. 
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